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We illustrate the main features of a recently proposed method based on ensemble density func- 
tional theory to divide rigorously a complex molecular system into its parts [M.H. Cohen and A. 
Wasserman, J. Phys. Chem. A 111, 2229 (2007)]. The illustrative system is an analog of the 
hydrogen molecule for which analytic expressions for the densities of the parts (hydrogen "atoms" ) 
are found along with the "reactivity potential" that enters the theory. While previous formulations 
of Chemical Reactivity Theory lead to zero, or undefined, values for the chemical hardness of the 
isolated parts, we demonstrate they can acquire a finite and positive hardness within the present 
formulation. 



1. INTRODUCTION 

In a series of recent papers [H, [H, [![, two of us have 
developed a rigorous method for dividing a complex sys- 
tem into its parts based on density-functional theory 
[H IS B 0, HI- The underlying theory, partition-theory 
(PT), was used to construct a formulation of chemical 
reactivity theory (CRT) [3] which, for the first time, is 
consistent with the underlying density-functional theory 
and is richer in structure than the preexisting CRT 




| , a sharp definition of the individual parts 
into which the whole system is partitioned is achieved 
first by selecting the nuclei of each putative part and 
maintaining these in the positions in which they occur in 
the whole and then requiring that the sum of the electron 
densities of the parts, each of which is treated as though 
isolated, add up exactly to the electron density of the 
whole (the density constraint). The electron densities of 
the parts are then to be determined by minimizing the 
sum of the density functionals of the individual parts with 
respect to the densities of the parts subject to the densit 
constraint. The density functional used, that of ref. 
(PPLB), allows for the existence of noninteger numbers of 
electrons on each part, necessary e.g. for the definitions 
of electronegativity (l2| and hardness [l3|, key indices 
of chemical reactivity [3j , and for incorporating covalent 
bonding between inequivalent parts. 

The minimization proceeds via a Legendre transforma- 
tion, which introduces a reactivity potential vr(t) as the 
Lagrange multiplier of the density constraint. Thus, the 
formalism can become computationally complex. First 
the electron density of the whole system must be deter- 
mined. Then, the densities of the parts must be deter- 
mined simultaneously with vr, all of which is required 



to set the stage for the determination of mutual reactivi- 
ties between parts, though certain self-reactivities can be 
determined for each species alone without reference to a 
larger system Q. 

Accordingly, in the present paper, we develope the par- 
tition theory in detail for an extremely simple system 
to exhibit its main features explicitly. The illustrative 
system is an analog of the hydrogen molecule in which 
the electrons move in one dimension along the molecular 
axis without interacting, and the nuclear Coulomb poten- 
tials are replaced by attractive delta-function potentials. 
As a consequence of these extreme simplifications, many 
quantities of interest can be determined analytically in a 
transparent manner, including the electron density of the 
molecule, of its parts (the "atoms"), and the reactivity 
potential at all internuclear separations. 

In Section 2, the model is defined and the molecular 
density obtained. In Section 3, the parts are defined, 
shown to have one electron each, and a polar representa- 
tion for their wave functions found which facilitates the 
minimization. In Section 4, the minimization is carried 
out, resulting in an Euler equation for the polar angle 
(3(x) of that representation. (3{x) is found in Section 5 
and used to determine the reactivity potential vr in Sec- 
tion 6. The principle of electronegativity equalization 
formulated in refs.Q and [|[ is shown to hold in Section 
7. Also in Section 7, the hardness [|| of the isolated H 
atom is calculated, shown to be nonzero, and correlated 
with the strength with which its electron is bound. Thus, 
despite the fact that the model is a caricature of the real 
system, meaningful features of the partition theory are 
indeed illustrated by it, as discussed in the concluding 
Section, 8. 
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2. 1D-H 2 ; INDEPENDENT ELECTRONS 
MOVING IN ATTRACTIVE 5-FUNCTION 
POTENTIALS IN ONE DIMENSION 

Our task is to partition an analog of the H2 molecule 
in which two electrons move independently in (5-function 
nuclear potentials in one dimension into parts, analogs 
of H atoms. Each H atom has, by symmetry only one 
electron, so the need for the PPLB density functional 
is avoided. Indeed no explicit use of density-functional 
theory is required for either the molecule or the atoms. 
The ground-state wave function ipo and energy Eg of an 
isolated H atom are (atomic units are used throughout): 



tpo{x) 



Z 2 /2 . 



(2.1) 
(2.2) 



In Eq. (|2.ip . (-Z) is the strength of the (5-function po- 
tential. To draw the analogy closer to real hydrogenic 
atoms, one could equate Z to the nuclear charge. 

The ground-state energy E{N = 1) of one electron 
moving independently in the two 5-function potentials 
centered at x = ±a is E(N — 1) = —k 2 /2, where k 
satisfies 



k = 2Z/(1 + tanh«a) 
The corresponding wavefunction is: 



= B C - 



cosh Ka 



\x\ > a 
\x\ < a 



where 

B = k 1 / 2 



1-1/2 



1 + 



cosh na 



tanh / 



(2.3) 



(2.4) 



(2.5) 



Note that k — > 2Z as a — ► (united atom limit) and 
k — ► Z as a — ► 00 (separated atom limit) . 

The two-electron molecular electron density is given 
by: 



um{x) = 2 \ip M (x)\ , 
and the total energy of the molecule is 

E M (N = 2) = 2E M (N = 1) = -k 2 



(2.6) 



(2.7) 



where N is the number of electrons in the molecule. The 
chemical potential of the molecule is therefore 

(J-m — E{2) — E(l) — E{\) — —k 2 /2 . 
3. PARITY DECOMPOSITION 



(2.8) 



We now partition the molecule into two parts a — 
1,2, each having a real one-electron wave function ip a , 
localized around —a and +a respectively, so that hm{x) 
is given by 



n M (x) = n%(x) + n 2 (x) , 



(3.1) 



where n a (x) is the electron density of each part a = 1, 2 
treated independently. The "atomic" wavefunctions are 
given by: 



■4> a (x) = \Jn a (x) . 
They are mirror images of each other, 
ip 2 (x) = ipi(-x) , 



(3.2) 



(3.3) 



and both are normalized. 

We now decompose the ip a into their symmetric, 
ip s (— x) — ips{x), and antisymmetric, ip a (~x) = —ipa(x), 
parts by a rotation within the function space they span, 



01 = A= (-0s + Ipa) , 1p2 = -J= {ips - 0a) 
0s = — (01 + 02) , 0a = — (01 - 02 N 



V2 



V2 



(3.4) 
(3.5) 



The rotation leaves "lengths" within the space invariant 
so that 



n M = ip 2 s +ipl 



(3.6) 



We next introduce j3 = (3(x), a polar angle in the function 
space, 



0s 



Inu cos (3 , ip a — y/n,M sin (3 



so that 



ipi,2 = y/n M /2 (cos j3 ± sin 0) 



(3.7) 



(3.8) 



Because the tp a are non-negative, |/3| cannot exceed 7r/4. 
Furthermore (3 must be an odd function of x, to ensure 
Q is also odd. This also guarantees normalization of ip a . 



4. THE EULER EQUATION FOR /3(x) 

To apply PT @,Q , begin with the original Hamiltonian 
H=\ /Zp^~ Z zZ [S(xi - a) + S(xi + a)] . (4.1) 



i=l,2 



i=l,2 



Then divide the system into overlapping regions, each 
with a given number of electrons. In this case, we choose 
one electron on the left, and the other on the right. Thus 
we have two 1-electron problems: 

2 

H a = ^+v a , v h2 = -ZS(xTa) . (4.2) 
The PT problem is to minimize 

£={ih,H 1 ih) + {ih,H2ih) , (4.3) 

subject to normalization of the wavefunctions, but also 
to the constraint that the total density equal the origi- 
nal molecular density, Eq. (|3.ip . (Without the latter con- 
straint, we'd obviously find ip\ t 2 = 0o(^ = T a ))- lu the 
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polar representation of Sec. 3, both density and normal- 
ization constraints are automatically satisfied, so the par- 
tition problem becomes simply minimizing £ as a func- 
tional of (3. That functional is 



dx 



1 



\^-\n' M + n M {0'f 
4 n M 2 



1 



+ [Oi + «2) + ("i - «a) si 112 / 3 ! f -( 4 - 4 ) 



Varying it yields 
S£ = I dx{n M 0'S0' + (v% - v 2 )n M cos 2/36/3} . (4.5) 



Integrating by parts, as usual, leads to 



5£ = 2n(3'5(3\ x x Z 



x— +oo 



+ {v 1 -v 2 )n M cos20}50 . (4.6) 



For £ to be stationary with respect to arbitrary varia- 
tions 8/3 of (3, both terms contributing to 5£ in Eq. (|4.6|) 
must vanish. The Euler equation which results from the 
vanishing of the second term in Eq. (|4.6|) is 



d/3 
dx 
dp 



dx V M dx 



d_ 

dx 



dx 



+ (vi -v 2 )n M cos 2/3 = , (4.7) 

+ Z (5(x - a) - S(x + a)) x 

xn M cos2/3 = . (4.8) 



The vanishing of the first term in Eq. (|4.6j) sets the bound- 
ary condition at infinity on the Euler equation H4.8[) . 
There are two possibilities, the vanishing of j3' at infinity 
or the fixing of (3 there so that 8/3 must vanish. As we 
shall see in Section 5, imposing the latter results in an 
unacceptable divergence in f3' at infinity. We therefore 
impose the boundary condition 



(3'{x) = , \x\ = oo . 

5. SOLVING FOR (3(x) 
Eq. (j4~8)) becomes 



~r ( nM ~r ) = ' \ x \ ^ a ) 

da; \ dx 



(4.9) 



(5.1) 



subject to the boundary conditions Eq. (|4.9[) and 

P(a~) = [3(a+) = f3 a 



0'(a-)-0'(a+) = Zcos2/3 a 

0(-a+) = P(-a-) = -0 a 
P'i-a') - P'{-a+) = Zcos2p a 



x = a , (5.2) 
x = —a . (5.3) 



The general solution of (|5.1[) is 



dP(x) 
dx 

p{x) 



Cl 



n M (x) 
dx' 



Cl 



n M (x) 



C2 



(5.4) 
(5.5) 



where ci and C2 are constants. As implied above in Sec- 
tion 4, if ci does not vanish 0' diverges exponentially at 
infinity, according to Eq. (|5.4|) , because um goes exponen- 
tially to zero, so, in accordance with Eq. (|4.9[) . c\ vanishes 
for |x| > a, and P(x) is constant there, 



P{x) — p a , x > a 
= -p a , x < -a 



For |x| < a we can rewrite Eq. (|5.6p as 

ci 



0(x) = 
which implies that 



dx' 



n M (x') 



0a 



a = - 2 



dx- 



ci 



n M (x) 

From (|5.8p we can relate ci to a via Eq. (|2.6[) , 

2nB 2 p a 



cosh na tanh kq 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



Inserting (|5.9p for ci into Eq. (|5.4[) and the result into the 
BC (|5.2p or (|5.3p produces an equation for p a , 
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Pa = — sinh 2na cos 2p a 
2k 



(5.10) 



Inserting Eqs. ()5.9p and 
markably simple result 



276]) into Eq.(f5~7l) yields the re- 



tanhra 

PW = 7 T Pa 

tann kq 



< \x\ < a 



(5.11) 



Eas.([5751). ([5710]) . and ([57IT]) . together with Eq. 
provide a complete analytic solution for /3(x) and through 
Eqs. (|3.5p and p.8p for the ip a - In FigfT] we show nM, 
ni and n 2 vs. x for Z — 1 and a = 1. We see that 
each localized density spreads into the neighboring re- 
gion, and looks quite similar to an atomic density. To 
see the differences from isolated atomic orbitals, in Fig [2] 
we make the distance smaller (a = 0.3), and show the 
right-side "atomic" orbital ipi(x) (solid line) and com- 
pare it with the pure exponential orbital ipo (x) of Eq |2.1l 
(dashed line). The orbital i/>i resembles ipo and tends to 
it for large a, but is distorted with respect to it for small 
a. Its maximum is still a cusp at x = a, but it also shows 
a second cusp at x = —a. Since k > Z always fEq. (|2.3p ). 
and either ipi or ■0 2 is proportional to ipM fo r 1^1 > a i 
where P = p a is constant, the PT atomic densities and 
orbitals decay more rapidly than isolated atoms. Since 
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FIG. 1: Molecular density um{x) (solid), and "atomic" den- 
sities ni(a;) and 712(0;) (dotted) for Z — 1 and a = 1. 




FIG. 2: Right-side "atomic" orbital ipi(x) (solid) and pure 
exponential orbital tpo(x) (dashed) for Z = 1 and a = 0.3. 



their normalization is the same, this in turn means en- 
hanced density between the 'nuclei', due to bonding. In 
FiglHl we show j3(x) for Z = 1, and a = 0.1, 1, and 10. 
Qualitatively, from Eq. (|5.11|) . 



(3{x) ~ f3 a — , x < min(l/K, a) 
a 



= Pa 



x > min(l/K, a) 



and if Za >> 1 (large separation), (3 a ~ 7r/4 while if 
Za << 1 (small separation), /3 a ~ a. The interpretation 
of these results is given in terms of (|3.8|) , outside the bond 
region. If /3 a is small, both 'atoms' share the density in 
each outside region. But if (3 a is close to 7r/4, each atom 
dominates on its own side, consuming the entire density 
there. 



1 

it/4 



-jc/4 
-1 













/ 0=10 
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a =0.1 













x/a 



6. THE REACTIVITY POTENTIAL 

The one-electron wave functions ipi(x) and ip2(x) are 
not eigenstates of the part-Hamiltonians Hi and H 2 of 
Eq. (|4.2|) . The natural question arises: What are they 
eigenstates of? The partition theory of refs. [l], 0, Q 
dictates that they are eigenstates of the modified single- 
electron Hamiltonians = p 2 /2 + V a , a = 1,2: 



Ipct = VMlpa , 
V Q = V a + Vr 



a = 1,2 



(6.1) 
(6.2) 



where the eigenvalue, regardless of the part a, is precisely 
equal to the molecular chemical potential fiM of Eq. (|2.8[) . 
The potential Vr{x) is the reactivity potential that we 
now construct explicitly. Summing over a and dividing 
by ipi + 1P2 yields a symmetric expression for vr, 



1 



vr = hm- 



V 



^1+^2 2 



(ipl + ^ 2 



Vflpl + V 2 l/J2 



(6.3) 



ipl and ip2 can be reexpressed in terms of ip s and ip a , 
Eq.([331). Noting that 



n M 



(6.4) 



using Eq. (|3.7p for tp Sia , and taking the S- function char- 
acter of v a into account results in 



1 



Vr 



Mm 



2ipM cos/3 dx 2 



(ipM COS (3) 



+u 2 )(l + tan/3 Q ) 



The molecular wave 
Schrodinger equation, 

l d 2 iP M 
2 dx 2 



function ip 



M 



(6.5) 

satisfies the 



+ (Vl + V 2 )lpM = MmV'M 



(6.6) 



which can be used to transform Eq. (l6.5p to 
2 



t'R 




dipM d/3 d 2 
dx dx dx 2 



tan /3 + i-f- 



d/3 
dx 



+ ^(vi+v 2 )(l 



tan p a ) 



(6.7) 



Using Eq. (|6.4jl . the Schrodinger-like equation for (3, 
Eq. (|4.8p , can be rewritten as 



2 d^p M d(3 



1 , 

+ 2 (U1 



-0 m efe tfrc 

tj 2 ) cos2/3 a 



dx 2 







(6.8) 



Multiplying Eq. (|6.8p by tan/3, invoking the oddness of 
[3 and the <5-functions in vi and U2, and subtracting the 
result from Eq. (|6.7p yields for vr 



FIG. 3: p{x) vs. x, as given by Eq. (|5.1ip . for fixed Z — 1 and 
3 different values of a. 



1 fd(3\ 1, 



(1 + cos2/3 a )tan/3 a ] . 

(6.9) 
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FIG. 4: Reactivity potential vr, Eg. (|6.10[) for fixed Z — 1 and 
3 different values of a: a = 0.1 (upper panel), a — 1 (middle) 
and a = 10 (bottom). The (^-functions at ±a are indicated by 
arrows. 



Inserting our previous result for /3(x), Eqs. (|5.6l) and 
(15. lip into (|6.9[) yields an explicit result for ur, 



vr 



I1M01 9(a-\x\) 
tanh 2 na cosh nx 

+ I{v 1 +v 2 )[l-sm2(3 a ] 



(6.10) 



where 6(y) — for y < 0, 1 for y > is the Heaviside step 
function. Eq. (|6. 10[) shows that vr(x) vanishes for \x\ > a, 
has attractive ^-functions at ±a whose weights increase 
monotonically from to \Z as -Za decreases from infinity 
to zero, and has an attractive inverse cosh 4 (a;) component 
for \x\ < a. For the united atom case, Za | 0, V\ + 
= i>2 + = 2v\ simply reproduces the molecular 
potential, and ip± = V>2 = ipM as they should. Figure 
2] displays us. x for fixed Z = 1 and representative 
values of a. The reactivity potential is almost flat for 
small separations, a wide well in between the two atoms 
for intermediate separations, and a narrow well that is far 
from both atoms at large separations. Figure [5] displays 
the weights of the (5-function components of vr divided 
by Z vs. a. 

As shown in ref.gj], the Kohn-Sham (KS) HOMO 
eigenvalue of each part must be identical to the chem- 
ical potential of the whole in the added presence of vr. 
In our simple example, the KS potential of a part re- 
duces to the nuclear 5-function potential of one H atom. 
Adding vr to the nuclear potential must therefore trans- 
form the HOMO energy E , Eq.flHJ), of the isolated 
atom to the more negative HOMO energy of the molecule 
E(N = 1) = —k 2 /2, which is its chemical potential 
(Eq. (|2.8[) ). vr must be attractive to do that, which it 
is, from Eqs. (|6.9p and (|6.10p . In our simple example, 
vr makes the delta function of the atom more negative, 
adds the attractive inverse cosh 4 potential between the 



FIG. 5: Weights of the 5-function components of vr divided 
by Z as a function of a for fixed Z = 1, from the second term 
of Eq. (|6.10p . The inset shows f3 a vs. a 
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FIG. 6: Energy as a function of a, in atomic units, for the 
two lowest-energy solutions of Eq. (|6.ip . Z = 1 for this plot. 



atoms, and adds an attractive ghost delta function at the 
position of the other atom to force the wave function to 
decay sufficiently rapidly outside the molecule. 

In the limit of infinite separation v\ +vr, reduces to Vi 
and V2 + vr reduces to v 2 , except for \x\ < a, where the 
attractive potential 



vr{x) 



1 



16 cosh Zx 



\x\ < a 



(6.11) 



persists. This potential has at least one additional bound 
state, but with binding energy less than |-Eo|- Thus it 
is unoccupied, and does not affect our results. The a- 
dependence of this state's energy is shown for fixed Z in 
FigJB] For very large separation between the atoms, it is 
localized at the center of the inverse cosh 4 (a:) component 
of vr, but it rapidly delocalizes for smaller separations. 
In particular, for Z = 1, it is highly delocalized when 
a <^ 1.4, where it vanishes into the continuum. 



7. SUSCEPTIBILITY AND HARDNESS 

Having found the reactivity potential, we now illus- 
trate the construction of reactivity indices. In the CRT 
of ref . Q , each part a is represented by an ensemble of 
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PPLB type containing contributions with only two in- 
teger electron numbers, p a and p a +i- The principle of 
electronegativity equalization is expressed as the equal- 
ity of the chemical potential of each part in the presence 
of the reactivity potential, n^, to the chemical potential 
of the molecule, [i m , 



Ma/ 



Va 



(7.1) 



The /i^ are defined as the difference between the ground 
state energies of a for p a + 1 and p a electrons in the 
presence of vr , 



and similarly for [1m 

Hm = Em(N m ) ~ E M (N M - 1) 



(7.2) 
(7.3) 

In our simple example, /im is given in Eq. (|2.8|) . The 
relevant value of p a is zero, so that /z^ is just E^(l), the 
lowest eigenvalue of 



Hi 



Vr 



(7.4) 



with H a given by Eq.([3~2"]l and v R by Ecu flTTD)) . The ex- 
plicit construction of vr in Section 6, not possible in gen- 
eral, guarantees that Eq. (|7.ip and therefore electroneg- 
ativity equalization holds. In the general case, a modi- 
fication of the Car-Parrinello scheme [13, HH guarantees 
electronegativity equalization. 

The susceptibility of part a measures the response of 
the density of part a to a small change in the potential 
V a of Eq.flBl]): 



X a (x,x') 



8n a {x) 



SV a (x>) ■ 

For 2 electrons, it is simple to show that 

X a (x 7 x') = -2Tp a (x)G a (n M ;x,x')ip a (x') 



(7.5) 



(7.6) 



where Q a {^M] x,x') is given by the E — > /im limit of: 



G a {E; x, x') = G a (E; x, x') 



ipajx^jx') 
E - Hm 

and G a is the Green's function for part a: 

-i -i 



G a (E;x,x') = 



e ->'t 



(x,x') 



(7.7) 



(7.8) 



Figure [7] shows the susceptibility of the right "atom" for 
various interatomic separations when the perturbing po- 
tential is added at xq — 3 (the numerical calculations 
were done as described in the Appendix) . Electrons flow 
away from xq, building up a peak at xq (positive because 
of the minus sign in the definition of \a, Eq. (|7.5|> ). and 
a negative peak at the closest maximum of the charge 
density, i.e. at a. With the analytic Green function of 



a = l 
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FIG. 7: Susceptibility x( x o,x) of the right "atom" obtained 
from Eqs, (17.6|l - (|7.8[) . as indicated in the Appendix, when xo 
is set to 3 a.u. Each panel corresponds to a different value 
of the internuclear distance, a. The lower-left panel shows 
\(xo,x) when a is just below (solid) and just above (dotted) 
x . 



an isolated "atom" [20| and Eqs.([71)|)-CLl]), Xa can be ob- 
tained analytically in the large-separation limit: 



Z(\x\ + \x'\) 

'} 



e -Z{\x\ + \x'\)\ e -Z\x'\ (j<^ 



We now construct the susceptibility of the whole sys- 
tem, \Ri by adding together the susceptibilities of the 
parts, 



Xr{x,x') = ^2xa(x,x') 



(7.10) 



The inverse of yr determines the hardness matrix r\ a {3 as 
shown in refs.Q] and Q: 



Vaf3 



dxdx' f a (x)x R 1 (x, x')fp(x') 



where the Fukui function of part a, f a (x), 

dn a (N a ,x) 



fa(x) 



dN a 



(7.11) 



(7.12) 



is simply equal to ipai x ) for 2 non- interacting electrons, 
since n a (N a: x) = Natp^x) (see also ref.[l6j]). Thus, we 
have 



Va/3 



dxdx'ipl(x)xR 1 (x,x')ipj(x') 



(7.13) 



Figure [8] shows the self-hardness rj aa for an isolated H- 
"atom" , as a function of Z. The constancy of the hard- 
ness for large Z can be understood qualitatively as fol- 
lows. The inverse susceptibility has units of energy times 
length squared. When Z is large, it establishes a length 
scale inversely proportional to Z, and an energy scale 
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proportional to Z 2 , so the Z-dependence cancels out in 
the inverse susceptibility. To obtain the hardness, we 
multiply Xr 1 on the left and right by the Fukui function, 
which has the dimension of inverse length. Integrating 
over position on the left and right then cancels out the 
Z-dependence arising from the Fukui functions, and the 
result is a Z-independent hardness. 



8. CONCLUSIONS 

Despite the extreme simplicity of the 1D-H2 model 
analyzed here - two non-interacting electrons moving 
in ID under the influcence of two equivalent attractive 
delta-function potentials - that model allows us to illus- 
trate the essential features of our partition theory and of 
key indices of our chemical reactivity via straightforward 
analysis and easy computations. 

We have shown that the electron density of the 
molecule can be decomposed exactly into a sum of atomic 
densities, a rigorous solution of the "atoms- in-molcculcs" 
problem [l7j ]. 

Electronegativity equalization [18| is built into the par- 
tition by the symmetry of the problem, so this homonu- 
clear model does not illustrate that principle as well as 
a heteronuclear model would. Nevertheless, the current 
example does illustrate a key feature of the new CRT, the 
chemical context dependence of the reactivity indices, 
in this case the electronegativity of a part, introduced 
through the presence of vr in the Schrodinger equation 
for ip a , cf. Eq. (|7.4p . It also demonstrates that the reac- 
tivity potential remains finite as two atoms separate, but 
has no effect on the partitioning after separation. 

Another serious shortcoming of the earlier formulations 
of DFT-based CRT is the vanishing of the hardness. We 
have shown explicitly here that the self-hardness, as de- 
fined in [H, of an isolated "atom" is positive. Interest- 
ingly, the hardness saturates as the ionization energy of 
the "atom" increases, raising the very interesting ques- 
tion of whether such a saturation of hardness with ion- 
ization energy exists in real systems. For this model, 
a strong positive correlation between hardness and ion- 
ization energy exists only over the limited range of Z 
between 0.4 and 0.7. 
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FIG. 8: Self-hardness vs. Z in the separated-atom limit 
(atomic units). 



where x < = inf (x, a;'), x > = sup(x, x') 

W[lp a ,L,1pa,R.] = i>a,h{E , x)lp' aR {E , x) 

-iP' a}L (E,x)iP a , R (E,x) , (A.2) 
and the orbitals ip a ,L and ip a .R are solutions of 

ip a ,L, R { E i x) = E^ L<R (E, x) (A.3) 



y + V Q (x) 



satisfying left and right-boundary conditions, respec- 
tively: 



\ip a , L {E,x)\ | , x i -oo 

\lpa,R( E ' X )\ I i x t °0 



(A.4) 
(A.5) 



The potential V a (x) of Eq. (|A~3l) is given by Eq. (|6.2p . with 
the reactivity potential vr(x) of Eo. (|6.10p . The compu- 
tations of tp a ,L,R{E, x) were carried out at E = hm±AE 
with AE chosen for numerical convenience, i.e. large 
enough so that sup x x ,\G a ([iM ± AE)\ does not become 
so large as to be inconvenient on the one hand, and small 
enough so that \ [G a (pM + AE) + G a (jiM ~ AE)] does 
not differ significantly from its limit at AE J, 0. We then 
calculated Q a of Eq. (|7.7|) as: 



KB is supported by NSF CHE-0355405. 



APPENDIX: Numerical calculation of the 
susceptibility 



Ga(fJ.M;x,x') = - [G a {nM + AE;x,x') + G a (fi M - AE\x,x') 

(A.6) 



We first obtained G a (E;x,x') according to the well- 
known prescription [lj|: 

G 9 {E^) = 2'ts^<^S^A , (A.l) 
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